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Abstract. Recently, Chen et al. [New J. Phys. 13 (2011) 083018] presented 
experimental results, accompanied by quantum-mechanical analysis, showing that 
the quantum interference behavior of Bell states could be simulated in a modified 
Mach-Zehnder interferometer whose inputs are pseudothermal light beams obtained 
by passing laser light through a rotating ground-glass diffuser. Their experiments and 
their theory presumed low-flux operation in which the simulated quantum interference 
is observed via photon-coincidence counting. We first show that the Chen et al. 
photon-coincidence counting experiments can be fully explained with semiclassical 
photodetection theory, in which light is taken to be a classical electromagnetic wave, 
and the discreteness of the electron charge leads to shot noise as the fundamental 
photodetection noise. We then use semiclassical photodetection theory to show that the 
same simulated quantum interference pattern can be observed in high-flux operation, 
when photocurrent cross-correlation is used instead of photon-coincidence counting. 



PACS numbers: 42.50.Ar, 42.50.Dv 
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1. Introduction 

The polarization Bell states of a pair of quantized spatiotemporal electromagnetic modes 
A and B are 

^ s \H U V) B ±\V)AH) Bt {l) 
v2 

and 

s \«)AH)B±\V)AV)B t (2) 
a/2 

where \H)k and |V)j<-, for K = A,B, denote single-photon states of horizontal and 
vertical polarization, respectively. These states form a maximally-entangled basis for 
the two-qubit Hilbert space of single-photon states for modes A and B. As such, they 
are extremely important for applications such as quantum teleportation p], quantum 
superdense coding [2J, and quantum key distribution [3], as well as their fundamental 
role in the Clauser-Horne-Shimony-Holt (CHSH) inequality [3]. The standard approach 
to generating these states is to post-select the biphoton output from spontaneous 
parametric downconversion (SPDC) [5J, and the standard approach to verify their 
entanglement behavior is via quantum-interference measurements [5]. Recently, Chen 
et al. [6] reported an experiment that mimicked the quantum-interference behavior 
seen with an SPDC entanglement source using two independent pseudothermal light 
beam obtained by passing laser light through a rotating ground-glass diffuser. Their 
experiments were carried out in the low-flux regime using photon-coincidence counting, 
and they provided a quantum-mechanical explanation that ascribed their observations 
to two-photon interference, just as is the case for SPDC light. 

It has long been known that the semiclassical theory of photodetection — in which 
light is treated as a classical electromagnetic wave and the fundamental photodetection 
noise is the shot noise arising from the discreteness of the electron charge — produces 
quantitatively identical predictions to those obtained from quantum photodetection 
theory when the illumination is in a classical state, i.e., a coherent-state or a statistical 
mixture of such states. See [7] for a detailed review of this topic. Except for any 
excess noise it may carry, laser light is coherent-state light. Moreover, its propagation 
through ground-glass diffusers, free space, optical fibers, and beam splitters are all linear 
transformations, for which classical-state inputs yield classical-state outputs. It follows 
that there must be an explanation for the experiments reported in [6] that relies on 
semiclassical photodetection theory, i.e., one that only needs classical electromagnetic 
waves. Our purpose in this paper is to present that explanation. Furthermore, 
although we will begin with a treatment that applies to low-flux operation using 
photon-coincidence counting, our approach readily extends to high-flux operation using 
photocurrent cross-correlation. 

The remainder of the paper is organized as follows. In Sec. 2 we describe the 
modified Mach-Zehnder interferometer — with photon-coincidence counting in the low- 
flux regime — that was employed in [6]. In Sec. 3 we introduce our classical-light model 
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for this interferometer, and use it to derive the singles and coincidence rates as functions 
of the interferometer's differential time delay and its polarization-analysis angles. Here 
we will show that our results explain the observations reported in [6]. Finally, in Sec. 4, 
we indicate how our theory can be extended to high-flux operation with photocurrent 
cross-correlation, and we provide some concluding remarks about the implications of 
our work. 

2. Modified Mach-Zehnder interferometer with pseudothermal Inputs 

The configuration for the experiment from [6] is shown in figure 1. A continuous- wave 
mode-locked Ti:sapphire laser operating at A = 780 nm wavelength with 78 MHz pulse- 
repetition frequency and a r p ~ 150 fs pulse duration illuminated an interference filter, 
to somewhat increase the pulse duration, followed by a rotating ground-glass diffuser, 
to render the light spatially incoherent. The diameter D « 5 mm output beam from the 
diffuser was divided by the 50-50 beam splitter BS1, with the resulting output beams 
propagating d ~ 200mm (from the diffuser) to collection planes, one of which could be 
offset, longitudinally, by 5. Each collection plane contained the tip of a single- mode fiber, 
whose transverse coordinates, p + and p_, satisfied \p + — p_\ ^> £ c , where £ c 31 /im 
is the correlation length of the speckles cast in these planes. The fibers routed the light 
they collected to polarizers PI and P2, set for orthogonal polarizations that we shall take 
to be x and y, respectively, before entering the 50-50 beam splitter BS2. The outputs 
from BS2 then underwent polarization analysis, by analyzers Al and A2 set for angles 
6a and 6b with respect to x, prior to single-photon detection. Computer processing 
completed the experiment by averaging the detector outputs over many pulses to obtain 
the singles rates Sa and Sb and the coincidence rate Cab- 

Chen et al. used quantum-mechanical analysis to show that Cab consisted of a 
background term plus a quantum interference term that is proportional to sin 2 (6a — 6b) 
when 5 = 0. Their background term arises from what in a related prior experiment 
[8] they called "self-intensity correlations." These correlations can be measured by 
summing the coincidence rates when a beam block is placed in front of the fiber tip at 
p + and when a beam block is placed in front of the fiber tip at p_. Subtracting the 
self-intensity correlations from the full coincidence rate obtained at 5 = then yields, 
according to theory, a unity-visibility quantum interference pattern. That background- 
subtracted quantum-interference pattern can then be used to obtain a violation of the 
CHSH inequality, in the same manner that is done [5] for the polarization-entangled 
signal and idler obtained from an SPDC source. The experimental data from j6] bears 
out this theory: the authors report 93.2 ±5.1% visibility in their background-subtracted 
quantum- interference pattern at 5 = 0. 

Chen et al. do not attempt to explain their experiments with a classical-light 
model, i.e., with semiclassical photodetection theory. In their previous work on a related 
anticorrelation experiment [Sj, they claimed there was no classical-light explanation for 
their results. We, however, have shown that such is not the case [9]. In particular, 
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Figure 1. Modified Mach-Zehnder interferometer with pseudothermal inputs, after 
[BJ. GG is a rotating ground-glass diffuser. IF is an interference filter. BSI and 
BS2 are 50-50 beam splitters, d is the nominal ground-glass to fiber-tip propagation 
distance, and S is a variable offset in one arm. PI and P2 are polarizers that select 
orthogonal polarizations that we shall take to be x and y, respectively. Al and A2 
are polarization analyzers set for angles 9 a and Ob, respectively. Da and Db are 
single-photon detectors 



we reported a straightforward classical-field analysis that reproduced the essential 
characteristics of the anticorrelation observed in [8]. Our demonstration is especially 
important for the following reason. Were Chen et al. [8] correct in asserting that their 
anticorrelation measurement could not be explained in this manner it would present 
quantum optics with a major conundrum: either laser light that has undergone linear 
transformation is not in a coherent state or a random mixture of coherent states, or the 
quantum and semiclassical theories of photodetection can make different quantitative 
predictions for the measurement statistics of classical-state light. In the present paper, 
we shall extend our scalar-wave model from [9] to provide a classical-light explanation 
for the simulated Bell-state reported in [6]. 

3. Singles and coincidence rates 

The linear velocity of the ground-glass diffuser where it was illuminated in [8] was 
~0.8m/s, so that for the r p ~ 345 fs and r p ~ 541 fs post-IF pulse durations considered 
therein it is fair to say that the ground glass was completely stationary while a single laser 
pulse propagated through it. We shall assume that to be the case for the experiments in 
[6]. The differential time delay 5t = S/c corresponding to the longitudinal variation S was 
±4 ps in [6] and we will assume that the photodetectors employed therein had the same 
T ~ 1 ns coincidence gate as in [8]. Hence, with r p < 1 ps, we have that \t p ± 8t\ <T. 
Because the fibers are single mode, and because their transverse coordinates have been 
set to ensure that the light beams they collect are uncorrelated, we will assume (cf. [9]) 
that E+(<) and E_(t), the positive-frequency classical fields emerging from polarizers 
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PI and P2, respectively, are given by 

E + (t)=v + f(t + St/2)e~^%, (3) 

and 

E-(t)=v-f(t-6t/2)e-**> t i v . (4) 

Here: v + and v _ are independent, identically distributed, zero-mean, isotropic, complex- 
valued Gaussian random variables with common mean-squared strength 

(\v + \ a ) = (\v.\ 2 ) = N, (5) 
representing constant-in-time speckle; 

is a transform-limited Gaussian pulse normalized to satisfy 

W(*)| 2 = l, (7) 

with r p being the post-IF pulse duration; and i x , i y are orthogonal unit vectors along 
the x and y directions. Note that in and (j3J) we have chosen ^/photons/s units, so 
that the average energy in either field is Nhuo, i-e., N is their average photon number. 

Because Chen et al. used photon-coincidence counting, their experiment was 
necessarily in the low-brightness regime wherein N ^ 1 prevails. As a result, we can say 
that the singles rates (counts/gate) and the coincidence rate (coincidences/gate) obey 

m 

fT/2 

S K = n dt(\E K (t)\ 2 }, iorK = A,B, (8) 

J-T/2 



and 



fT/2 nT/2 

C AB = V 2 dt du(\E A (t)\ 2 \E B (u)\ 2 ). (9) 

J-T/2 J-T/2 

In these expressions: 77 is the photodetectors' quantum efficiency; the complex envelopes 
of the fields that illuminate detectors D A and D B are 

cos(6 A )v + f(t + St/2) + sm(9 A )v.f(t - 5t/2) 
E A (t) = -j= , (lUj 

and 

cos(6 B )v + f(t + g/2) - sm(e B )v- f(t - 5t/2) 

&B{t) = , (11) 

where we have suppressed all propagation delays except for the interferometer's 
differential delay St. Using the statistics of v + and v _, along with the disparities between 
t p , St, and T, it is easy to evaluate the singles rates and the coincidence rate. 
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The statistical independence of the v + and t>_, and their mean-squared strengths, 
immediately gives us 

n ( r T/2 

I I /L. |2\ \ / j- I ft, , 



^ = 7^ (Kr)cos^(^) / dt\f(t + 5t/2)\ 

- \ J-T/2 

+ (K| 2 )sin 2 (^) / ^|/(t-5t/2)| 2 (12) 

J-T/2 J 

for AT = A, B, (13) 

where the approximation follows from |r p ± 5i/2| < T. Similarly, for the coincidence 
rate, the statistical independence of v + and t>_ leads to 
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Cab = -7- 1 44(|^| 4 ) dt |/(t + )| 2 du \f{u + )\ 



T/2 J-T/2 



rT/2 r T/2 

+ 44<KI 2 )<M 2 > / dt\f(t + )f / du\f(u_)\ 2 

J-T/2 J-T/2 
rT/2 r T/2 

+ 44(KI 2 XM 2 > / dt\f(t„)\ 2 / ^|/K)| 2 

J-T/2 J-T/2 
rT/2 rT/2 

-2c A c B s A s B (\v + \ 2 )(\v4 2 ) / dt duRe[f*(t + )f*(u-)f(t-.)f(u + )] 

J-T/2 J-T/2 
rT/2 r T/2 \ 

+ 44(K! 4 ) / dt\f(t^)\ 2 du\f(u„)\ 2 ), (14) 

J-T/2 J-T/2 J 

where ck = cos(6 l x) and sk = sm^Ox), for K = A, B, w± = w ± 5t/2, for w = t,u, and 
the first and last terms on the right are the self-intensity correlations. Following [6], 
we shall suppress these self-intensity correlations — they can be found by first measuring 
the coincidence rate with E + (t) blocked and then measuring the coincidence rate with 
E_(t) blocked — by subtracting them from Cab and focus our attention on 

2^2 / rT/2 rT/2 

Cab = 1 —[c 2 a s 2 b / dt\f(t + )\ 2 du\f(u_)\ 2 

4 V J-T/2 J-T/2 

rT/2 r T/2 

+ c 2 B s\ dt\f(t.)\ 2 du\f(u + )\ 2 

J-T/2 J-T/2 

rT/2 rT/2 \ 

- 2cac b s a sb / dt duRe[f*(t + )f*(u-)f(t-)f(u+)]} , (15) 

J-T/2 J-T/2 J 

where we have used the mean-squared values of v + and t>_. Now, by exploiting 
\t p ± 8t/2\ <C T and some trigonometric identities, we can reduce the preceding 
expression to 

7] 2 N 2 



C AB ^ 



_ Qb) + sin(2g A )sin(2g fl ) (1 _ 



(16) 
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At this point we are ready to compare our classical-light theory with the 
experimental results from [BJ. There it is shown that 9 a = 9b = vr/4 leads to an 
ant icor relation dip of width ~r p in Cab as St is scanned through zero. Moreover, this 
anticorrelation dip has near- unity visibility, viz., 

m&x(C AB ) - mm(C AB ) 



max(C AB ) + nrin(CAB) 
From ffTBl) we find 



1. (17) 



AT 2 

C AB = V -^-(l-e^), (18) 

in agreement with those observations. Likewise, Chen et al. find that when 9a = 
and 9 B = 37r/4 there is a near-unity visibility correlation peak of width ~r p in Cab as 
St is scanned through zero. From (fT6l) we get 

C AB = ^f(l + e-^), (19) 

for this case, as was found in the experiment. Finally, at St = 0, [BJ reports simulated 
Bell-state quantum- interference when 9a — 9b is varied. That characteristic is present 
in (fT6|) . where we find 

C A B = 1 ^sm 2 (9 A -9 B ) (20) 

when St = 0. Collectively, these matches between the Chen et al. experimental 
results and the predictions of our classical-light theory are consistent with their 
experimental configuration's illuminating its detectors with classical-state light and the 
quantitative agreement of the quantum and semiclassical theories of photodetection 
for such illumination. Because light is quantum mechanical, and photodetection is 
a quantum measurement, the quantum explanation of their work is of course the 
fundamental one, but we have just shown that it is not necessary to invoke a quantum 
description of the light to explain their results. 



4. Discussion 



Our work in Sec. 3 presumed low-flux (single-photon) operation, specifically iV <C 1, in 
deriving the singles and coincidence rates for the figure 1 experiment. Chen et al. raised 
the question of what would happen if the illumination in this experiment was performed 
in the high- flux regime, wherein N ^> 1 prevails [6j. With our classical-light theory, 
it is easy to answer that question. Suppose that E-t(t) are as given in Sec. 2 except 
that N ^> 1. Also assume that the single-photon detectors in figure 1 are replaced by 
shot-noise limited PIN photodiodes, whose photocurrent outputs are lAif) and is(t), 
respectively, and that the singles and coincidence-rate measurements are replaced by 
average photocurrent and photocurrent cross-correlation measurements, = (^(0)), 
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for K = A,B, and lf B = (i A (0)i B (0)). We then have that [7] 

JW = TL ((K| 2 > cos 2 (^) | rfr \f(r + 5t/2)| 2 />(-r) 

+ <|w_| a >sm a (te) jdr\f{T-5t/2)\ 2 h{-T)^ , (21) 

where g is the electron charge and h(t) is the photodiodes' baseband impulse response, 
which is normalized to satisfy jdth(t) = 1. To make the connection to the case of 
photon-coincidence counting, we shall ignore causality and take h(t) to be 

f 1/T, for \t\ < T/2 
hit) = { (22) 
I 0, otherwise, 

with T ~ 1 ns, so that we again have \r p ±5t/2\ < T. We then can reduce our expression 
for the average photocurrents to the following simple form: 

tf'*^. for K = A,B. (23) 
For the photocurrent cross-correlation we have 

™ a 2 n 2 ( f T/2 f T/2 

^1 = ^ 44(K| 4 ) / dt\f(t + )\ 2 du\f( U+ )\ 2 
\ *J — T l*2i J — T/2 

f T/2 rT/2 

+ 44<K| 2 )<K| 2 > / dt \f(t + )\ 2 / du \f{u-)\ 2 

J-T/2 J-T/2 
r T/2 pT/2 

+ 44(KI 2 XM 2 > / dt\f(t_)\ 2 / rf M |/K)| 2 

J-T/2 J-T/2 
rT/2 r T/2 

-2c A c B s A s B (\v + \ 2 )(\v4 2 ) / dt duRe[r(t + )r(u-)f(t-)f(u + )} 

J-T/2 J-T/2 
rT/2 r T/2 \ 

+ 44(K| 4 }/ dt\f(t.)\ 2 du\f{uJ)\ 2 ), (24) 

J-T/2 J 

where we have used (1221) for /i(t), and the first and last terms on the right come from self- 

(2) 

intensity correlations. Subtracting the self- intensity correlations from I AB , and using 



(2) iu , u . ium( . 

2„2 7\r2 / /-T/2 /-T/2 



to denote the result, we arrive at 



„2„2a/-2 / ri/v rl/A 

if B = q -^- 44 / ^l/(* + )l 2 / ^i/Mi 2 

\ J — -T/2 */ — T/2 



fT/2 /-T/2 

+ 44/ dt\f{U)\ 2 du\f(u + )\ 2 

T/2 J-T/2 
T/2 fT/2 



- 2c A c B s A s B dt duRe[f*(t + )r(u-)f(t-)f(u + )]) (25) 

J-T/2 J-T/2 / 

q 2 r] 2 N 2 



4T 2 



T/2 J-T/2 

2 _ ^ + sin(2^) sin(2^) _ 



(26) 
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Thus the same simulated Bell-state behavior seen in photon-coincidence counting in 
low-flux operation will be present in the photocurrent cross-correlation in high-flux 
operation. 

Let us conclude with a brief meta-lesson that can be gleaned from what we 
have done. Chen et al. [6] use quantum analysis to show that photon-coincidence 
counting in the figure 1 configuration — after subtraction of self-intensity correlations — 
mimics Bell-state quantum interference. Because classical-state light is employed 
in their experiments, their observations have an equivalent quantitative explanation 
from semiclassical photodetection theory. Moreover, as we have just seen, that same 
mimicking of Bell-state quantum interference will be present in photocurrent cross- 
correlations, taken in high-flux operation, after self-intensity correlations are subtracted. 
In the high-flux regime, it is clear that semiclassical theory ascribes the interference 
pattern to intensity-fluctuation correlations between |£U(i)| 2 and \EB(t)\ 2 . Inasmuch 
as the experimental results from |6j cannot distinguish between the quantum and 
semiclassical theories, it is disingenuous to claim that the figure 1 configuration generates 
Bell states, as opposed to simulating the quantum-interference pattern produced by 
detection of a Bell state. After all, Bell states are entangled states, and entangled states 
are nonclassical. Hence their photodetection statistics cannot be fully and properly 
quantified by semiclassical photodetection theory. In this regard we note that post- 
selection, which is used to identify Bell-state photon pairs produced by SPDC, is very 
different from subtracting the self-intensity correlations in the figure 1 experiment to 
see Bell-state quantum interference. This is because post-selection is just selecting the 
occurrence of a photon coincidence, i.e., it can be performed on a pulse-by-pulse basis, 
whereas subtracting the self-intensity correlations in figure 1 requires collecting averages 
of coincidences with both beams present at BS2 and with only one beam or the other 
present on that beam splitter. 
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